Tip Splittings and Phase Transitions in the Dielectric Breakdown Model: Mapping to 

the DLA Model. 
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We show that the fractal growth described by the dielectric breakdown model exhibits a phase 
transition in the multifractal spectrum of the growth measure. The transition takes place because the 
tip-splitting of branches forms a fixed angle. This angle is r\ dependent but it can be rescaled onto an 
"effectively" universal angle of the DLA branching process. We derive an analytic rescaling relation 
which is in agreement with numerical simulations. The dimension of the clusters decreases linearly 
with the angle and the growth becomes non-fractal at an angle close to 74° (which corresponds to 
?/ = 4.0 ±0.3). 



Fractal growth and patterns are common phenomena 
of nature. The prototype model for a mathematical de- 
scription of fractal growth is the diffusion-limited aggre- 
gation (DLA) model Q. The growth in this model is 
determined by the electric field (called the harmonic mea- 
sure) around the emerging fractal cluster. The harmonic 
measure possesses multifractal scaling properties and re- 
cently, a new insight into the behavior of this measure 
was presented by applying the method of infinitely con- 
voluted conformal mappings Q . In particular it has been 
demonstrated that there exists a singular behavior in the 
multifractal spectrum and this singularity is a signature 
of a phase transition in thermodynamics formalism of 
the spectrum 0. The phase transition in the DLA clus- 
ter occurs at a specific moment q c of probabilities of the 
growth. It has further been demonstrated that the geo- 
metrical properties of the DLA cluster, which gives rise 
to this transition, is the existence of a specific branching 
angle for each new offspring on this cluster |5| . For DLA 
this critical branching angle was found to be around 27°. 
Another approach which defines a characteristic angle in 
the DLA process, is to consider the stability of a finger 
growing in a wedge || . 

The DLA model has been nicely generalized to the di- 
electric breakdown model (DBM) || where the growth 
probabilities at a specific site of the cluster is determined 
by the value of the electric field (or the harmonic mea- 
sure) raised to a power ry, i.e. the growth measure at the 
interface follows, p v (s)ds ~ \E(s)\ v ds where DLA cor- 
responds to the case r\ = 1. With varying values of ij 
clusters of different geometry are grown each with their 
own characteristic properties of the multifractal spec- 
trum. These properties have not been outlined before 
and it is the purpose of this letter to examine in details 
possible critical points and phase transitions in the ther- 
modynamic formalism of the harmonic measure for the 
dielectric breakdown model, at varying values of 77. 

The structure of the clusters of the DBM model 
emerges from an on-going proliferation and screening, 
hence stagnation, of branches. In particular will the pro- 
truding branches create fjords in which the harmonic 
measure decreases rapidly compared to what happens 
around the tips. The multifractal properties are best 
studied using the recently proposed model of iterated 



conformal maps Q , since the deep fjords, numerically, are 
invisible to the original approach where random walkers 
are used as probes, see ||. The model is based on com- 
positions of simple conformal maps 4>\. e which take the 
exterior of the unit circle to its exterior, except for a lit- 
tle bump at e ie of linear size proportional to We 
shall here use the mapping introduced in Q which pro- 
duces two square root singularities which we refer to as 
the branch cuts, and the tip of the bump which we refer 
to as the micro tip. The composition of these mappings 
is analog to the aggregation of random walkers in the 
off-lattice DLA model. The dynamics is given by 



(1) 



The size of the n'th bump is controlled by the param- 
eter X n and in order to achieve particles of fixed size we 
have that, to first order, 
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The growth probability pi(s) at the interface of a DLA 
cluster of size n is, in the electrostatic picture, propor- 
tional to the electric field 



ds 



Pl (s)ds ~ \E(s)\ds ~ — . 



(3) 



In this case the measure on the unit circle is uniform. 
When the electric field is raised to the power 77, the mea- 
sure is no longer uniform, 



Pv (6)d6 ~ p v (s(6)) 



ds 
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Numerically we use the Monte Carlo technique intro- 
duced in II in order to choose 9 according to the dis- 
tribution Pri . We vary the number of iterations T such 
that for r) = 1.25, T — 50 and for T) = 4, T = 400. 

Below we consider both the growth measure and the 
harmonic measure. The growth measure is used to derive 
the multifractal properties whereas the harmonic mea- 
sure is used to determine the physical properties. Let 
us emphasize that during the growth we always use the 
growth measure. First, we consider the growth measure 
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and the phase transition in the corresponding multifrac- 
tal spectrum. 

The moments of the growth probability scales with 
characteristic exponents, the generalized dimensions, 

J p q {s)ds ~ (l/JQfo- 1 )^ ~ n (i-i)D q /D j (5) 

where n is the number of particles and R the linear size 
of the cluster. 

Numerically we approximate the integral on the left 
hand side by the sum of the field evaluated along the mi- 
cro tips of the bumps produced by the bump mappings. 
The field in DLA will for clusters of size 20000 assume 
values below 10 -20 and it follows from ^) that it is im- 
possible with the numerical precision on the unit circle 
(A9 ~ 10~ 16 ) to maintain a reasonable resolution in the 
physical space As = \E\A9 ~ 10 4 . We therefore use the 
resolution increasing approach introduced in [0] where 
one keeps track on the dynamics of the branch-cuts. 

An easy way to see the existence of a phase transition 
in the multifractal spectrum is to look at the distribution 
of p v sampled along the tips of the bumps. The distribu- 
tion will for the smallest values (below some cutoff value 
c) of prj < c scale with a characteristic exponent ^73^ ■ 
The value of j3 is calculated by reordering the N com- 
puted values of p v in ascending order. In other words, 
we write them as a sequence {p v (i)}i£i where / is an 
ordering of the indices such that p v (i) < P-qii) if i < j- 
We treat the discrete index i/N as a continuous index 
< x < 1 and therefore p v as a non-decreasing function 
of x, 

Pv = f( x ) ■ ( 6 ) 

Numerically we find that the function f(x) obeys a power 
law with an exponent (3 for values of x <C 1. From f(x) 
we calculate the distribution of p(p n ) by the transforma- 
tion formula 

P(P V ) ~ J 5{Pn f{x))dx = lr{x \ pv))l ~ p/ ■ (7) 

With the use of this distribution the moment integral (|5|) 
is rewritten as 

/>L />oo 
/ P q n ds = I p q v p{pr,)dpr, 
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= k J P^Pv' 3 dPv + J p q v p(p v )dp v , (8) 

where k is some normalization constant. The left integral 
in the final expression diverges whenever 

3 < <Zc - ~ (9) 

The phase transition in the multifractal spectrum of the 
growth measure takes place for the value q = q c . We shall 




FIG. 1: Part of a cluster grown with r\ = 3. The wedge 
structure at the bottom of the fjords is clearly seen and the 
opening angle observed along the aggregate is close to the 
angle predicted in ( |l3[ ) and shown on the figure. 

now argue that the value of /3, and therefore q Cl is inde- 
pendent of the value of 77. In other words the phase tran- 
sition in the multifractal spectrum of the growth measure 
is universal. 

In Q it was shown that the angle defining the branch 
splitting deep inside the fjords of DLA was given by a 
characteristic angle 7 C (1) = 7 C (?7 = 1). This characteris- 
tic angle was also identified as the reason for the phase 
transition in the multifractal spectrum of the harmonic 
measure of DLA. The electric field along the branches in 
a wedge with opening angle 7 scales like 

pi(x) oc \E(x)\ - ir 1 (10) 

and therefore in DBM the growth probability inside this 
wedge is given by 

p v (x) x\E(x)\^ x^~ 1} . (11) 

When the exponent r\ is introduced, the growth proba- 
bility (|ll]) inside the wedge is changed and therefore we 
see a corresponding change in the geometry, such that 
e.g. when r\ > 1 the field inside a wedge with opening 
angle 7 C (1) is effectively similar to that of a wedge with 
a smaller opening angle. 

The wedge structure at the bottom of the fjords of 
DLA is not affected by a change in the value of 77 see 
Fig. [l], the physical angle, however, is changed such that 
the effective angle remains fixed and equal to the angle 
observed in DLA. Generally we therefore have for the 
angle observed physically that j c (ri > 1) > 7 C (1) and 
7 C (77 < 1) < 7c(l)- The angle for a given value of 77 is 
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FIG. 2: The critical angle as function of r\. The line represents 
the critical angle obtained from the analytical expression (|l3|). 
Each dot represent the numerical average of 4-15 clusters of 
size n — 18000. The standard deviation of the individual 
points assumes values in the range between 1.5° and 3.4°. 



found by comparing the exponents in ( pd| ) and ( |To| ) 
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Fig. g shows together with numerical predictions (see 
below) how the critical angle varies with rj. The value 
of 7c(l) in (|l^) and used in the figure is determined nu- 
merically from 15 DLA clusters of size n=20000, 7c(l) = 
27° ±3°. 

The distribution ofp,, inside a wedge with opening 
angle 7 follows from (Ull), with a — tt/j, 



p(Pv) 



1 



( x (Pn)) 



rj(a-l)-l 



Pv 



(14) 



and if we compare the exponent with that of (]?]) and 
insert the critical angle from ([[2]) we find that 



7c(l) 



- 1 



(15) 



Therefore (3 in expression (|9j) is independent of 77. 

One way, numerically, to calculate the critical angles 
shown in Fig. His first to locate the regions where the 
distribution in (m) scales and afterwards perform a direct 
measurement in these regions. Such measurements are 
most likely rather inaccurate and therefore we turn to 



FIG. 3: The estimated dimension plotted versus the critical 
angle j c , see (fiM- The range of the critical angle corresponds 
to values of 77~Detween 0.75 and 3.5. The linear fit of the 
points intersects the line D=l at an angle y c ~ 74°. 



(JlO|). During the growth we apply the exponent of 77 as 
usual but once a cluster is grown we consider the har- 
monic measure p\ only. Similarly, as above, we find that 
inside a wedge 



p(Pi) ~ Pi 



(16) 



The exponent is compared to the one we calculate numer- 
ically from (^) and from this we find the critical angle as 
function of (3. Note that when we consider the harmonic 
measure and not the growth measure, p is a nontrivial 
function of 77 and the critical angle is given by 



7c(»7) 



(17) 



Another interesting observation is that the dimension 
seems to depend linearly on the critical angle as shown 
in Fig. ||. Note that the linear fit intersects the line of 
value equal one before the angle reaches 180 degrees, and 
therefore the growth is only fractal for rj below a finite 
value. 

The point of intersection is found at the angle -f c — 74° 
and using (|l^) we find that 

77 = 4.0 ±0.3, (18) 

in agreement with the results obtained in |Q and M . To 
sum up on these results, we rewrite the relation between 
the dimension and the critical angle in terms of the angles 
72 and 71 at which the growth becomes two- and one- 
dimensional respectively 



£>(7c) = 1 + 



7i ~7c 
71 - 72 



72 < 7c < 7i • 



(19) 
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The dimension can also by (13) be written in terms of ij 
and in this case the dependence will no longer be linear 
but be of a form similar to that observed in j7j. Due to 
the finite size of the bumps used in the growth we do 
not observe, as rj — ► + , that the critical angle of the 
fjords vanishes. The bumps will fill up the fjords, and 
the growth become two dimensional for a non-vanishing 
value of 77 and 72 ~ 5°. The fill up problem is also the 
reason why we have not been able to present data points 
for values of 77 below .75, because the deep fjords are filled 
up. 

In conclusion, we have studied the critical properties of 
the growth of clusters in the dielectric breakdown model. 
In particular, we have focused on the branching process 
and have found that the branching, on the average, oc- 



curs at a fixed angle which depends on the value of the 
characteristic parameter rj of the DBM model. The size 
of the angle in turn determines the phase transition point 
of the growth measure. We have derived an analytic ex- 
pression for the variation of the branching angle with rj 
and found excellent agreement with numerical data. Fur- 
ther, we have found a linear dependence of the dimension 
of the cluster with the value of critical angle. This linear 
dependence results in a prediction of the branching angle 
at the point where the growth becomes one-dimensional. 
It is found to be j c » 74° corresponding to r\ c » 4 in 
agreement with a result obtained in and ||. We are 
in the process of developing a scaling theory for the linear 
dependence of the dimension versus angle. 
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